Decoherence induced by a dynamic spin environment (II): 
Disentanglement by local system-environment interactions 
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This article studies the decoherence induced on a system of two qubits by local interactions with 
a spin chain with nontrivial internal dynamics (governed by an XY Hamiltonian) . Special attention 
is payed to the transition between two limits: one in which both qubits interact with the same site 
of the chain and another one where they interact with distant sites. The two cases exhibit different 
behaviours in the weak and strong coupling regimes: when the coupling is weak it is found that 
decoherence tends to decrease with distance, while for strong coupling the result is the opposite. 
Also, in the weak coupling case, the long distance limit is rapidly reached, while for strong coupling 
there is clear evidence of an expected effect: environment-induced interactions between the qubits of 
the system. A consequence of this is the appearance of quasiperiodic events that can be interpreted 
as "sudden deaths" and "sudden revivals" of the entanglement between the qubits, with a time scale 
related to the distance between them. 

PACS numbers: 03.65.Yz, 03.67.Lx, 03.67.Mn 



I. INTRODUCTION 

Understanding the process of decoherence [J, 0, 0| is 
not only important from a fundamental point of view 
but also essential to design good error correction strate- 
gies to prevent the collapse of quantum computers [4J. 
When a composite system (for example, a quantum com- 
puter with many qubits) interacts with an environment, 
the process of decoherence generates disentanglement be- 
tween the components of the system. Different qubits of 
a quantum computer would interact with environments 
that may sometimes be effectively independent, but in 
other cases the correlations in the environment may not 
be negligible. In such a context, it is important to under- 
stand the properties of the decoherence process and the 
signs in the decay of quantum coherence of an environ- 
ment which is not only dynamic (i.e., it has a nontrivial 
evolution of its own) but also spatially correlated. In this 
paper we make a step in this direction. 

In a previous paper we analyzed the decoherence in- 
duced on a single qubit by the interaction with an envi- 
ronment formed by a spin chain with an XY Hamiltonian 
Q . The way in which quantum coherence decays in time 
was found to depend strongly on the nature of the envi- 
ronment's Hamiltonian (see [5[ for a list of references of 
previous studies on spin baths). Here, we will consider 
a system formed by two qubits interacting with an envi- 
ronment modelled in the same way. Each qubit in our 
system will be put in local contact with separate sites of 
the chain. This situation is physically more realistic and 
interesting than the more familiar one of central qubits 
with homogeneous couplings to all the chain. We note, 
however, that our system is equivalent to a one-qubit 
system with non-local couplings. The decoherence pro- 
cess will induce a decay of the entanglement between the 
qubits that will depend on the distance between them. 
We will study different scenarios ranging from the case 
where both qubits interact with the same site, to the long 



distance situation which is similar to the action of two 
independent environments [6| . We will specially consider 
two limiting cases of very weak and very strong system- 
environment couplings, and see that the decoherence ef- 
fects differ substantially. In the strong coupling regime, 
the chain provides a medium through which the qubits 
can interact, so that the entanglement between them ex- 
hibits quasiperiodic events of "sudden deaths" and "sud- 
den revivals" [7[ with a time scale that depends on the 
distance between qubits. 

The structure of the paper is as follows: in Section 
HT1 we introduce the model, defining the Hamiltonians 
for the system, the environment, and the coupling be- 
tween them. We also present the main formulas we will 
use to determine the decay of quantum coherence. In 
Sections IIIII and IIVI we study the decay of the entan- 
glement between the qubits as a consequence of their 
non-homogeneous interactions with the environment, in 
the cases of weak and strong couplings with the chain. 
Finally, in Section [V] we summarize our results. 



II. THE MODEL 

A. The system, the environment, and their 
evolution 

We will study the decoherence induced on a system of 
two spin 1 /2 particles (which we shall call "the system" or 
"the qubits" ) by the coupling to an environment formed 
by a chain of N spin 1/2 particles. We neglect the self- 
Hamiltonian of the system, and take the Hamiltonian of 
the environment chain as: 
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where Xj , Yj , Zj denote the three Pauli operators acting 
on the j-th site of the chain, and periodic boundary con- 
ditions are imposed. The parameter 7 determines the 
anisotropy in the x — y plane and A gives a magnetic field 
in z direction (7 = 1 corresponds to the Ising chain with 
transverse field). This model is critical for 7 = 0, |A| < 1 
and for A = ±1, which corresponds (in the limit N — > 00) 
to a quantum phase transition from a ferromagnetic to 
a paramagnetic phase. Throughout the paper, we shall 
discuss the effects of this phase transition on the deco- 
herence of the system. 

Each of the two qubits is made to interact locally with 
a spin of the chain: the first qubit (5*1) with the spin 1 
(this choice is arbitrary because of the periodic boundary 
conditions), and the second qubit (S2) with the spin n. 
The interaction Hamiltonian is chosen as: 

H int = -9 (|l)(lk^i + |l)(l|s 2 ^) (2) 

with |0) and |1) the two eigenstates of the Z Pauli oper- 
ator. Thus, if the system is in state \ab) (a,b E {0, 1}), 
the environment evolves with an effective Hamiltonian 
H ab given by: 

H ab = H c - g{aZ x + bZ n ). (3) 

In this way, the effect of the coupling is a change in the 
effective external field of spins 1 and n, conditioned on 
the state of the system. 

We assume the initial state of the universe formed by 
the system and the environment to be pure and separable 
(i.e., the environment is not correlated with the system): 

Pse(0) = \ip)(tp\ <g> \E )(Eo\ (4) 

Our goal is to study the evolution of the reduced den- 
sity matrix of the two-qubit system (obtained from the 
state of the universe by tracing out the environment). 
Because of the special form of the Hamiltonian, the tem- 
poral dependence of the reduced density matrix p can be 
formally obtained as follows. In the computational basis 
of the two-qubit system (formed by the eigenstates of the 
Z Pauli operator of each qubit), p can be written as: 

p{t)= J2 Pab,cd(t)\ab)(cd\ (5) 

abed— 0,1 

The evolution of the elements of p is given by: 

PabMt) = PabMO){ E o\e iHBdt e~ iH, * t \Eo) (6) 

(we are using units such that h = 1). As the total Hamil- 
tonian commutes with Zs 1 and Zg 2 , the diagonal terms 
Pab.ab remain constant. On the other hand, each off- 
diagonal term in the reduced density matrix is modified 
by a factor smaller than 1. This factor is associated with 
the overlap between two states of the environment that 
correspond to the two different evolutions of the spin 
chain according to the different system states. To ana- 
lyze the decoherence induced by the coupling to the spin 



chain we will consider the square modulus of this factor, 
which (after Q) we shall call the Loschmidt echo: 

LabMt) = \(E \e lH ^e-* H ^\E )\ 2 . (7) 

This echo is simplified if we assume the environment to 
be initially in the ground state of one of the effective 
Hamiltonians involved. In such a case, one of the evo- 
lution operators in the expression acts trivially, and the 
echo is then equal to the survival probability of the initial 
state after being evolved with the other operator. 

To find the decay of coherence with time for this prob- 
lem we note that the Hamiltonians H a b of the chain 
can be mapped onto a fermion system by means of the 
Jordan- Wigner transformation |8j: 




Zj = 2c] Cj -l. (10) 



Using this, up to a correction term associated to bound- 
ary effects, the Hamiltonians can be written as: 

H ab = ~ yj(cjc J+ i + c] +1 c 3 ) + 7(4 c I+i + c J+^ c o) + 
3 

+Xf\2c]c 3 -l) (11) 

with X^ ab ^ = A + g(aSj^i + b5j :n ) (the extension to the 
case where the qubits interact with more than one site is 
trivial) . 

The Hamiltonians H ao depend quadratically on the an- 
nihilation and creation operators. Therefore they can 
be diagonalized by linear (Bogoliubov) transformations 
defining new creation and annihilation operators which 
we will denote as r)( ab \ rj^ ab \ Furthermore, as all these 
transformations are linear, the operators corresponding 
to different values of the labels (ab) can also be connected 
by Bogoliubov transformations. 

The factors modulating the off-diagonal terms of the 
reduced density matrix can be written in terms of these 
transformations. If the initial state of the environment 
\E ) is the ground state of the unperturbed chain Hamil- 
tonian _ffoo we can show (see Appendix 1) that: 

LabMt) = K^ole-^'lSo)! 2 = |det(g + he iAt )| 2 (12) 

Here A is a diagonal NxN matrix containing the energies 
of the normal modes of the Hamiltonian H a \>. g and h 
correspond to the transformation connecting the particles 
that diagonalize the unperturbed Hamiltonian _ffoo an d 
the effective Hamiltonian H ab : 

4 00) -Egrf b) +W (ab) (13) 

k 



In this way we can express the echo as the determinant of 
anJVxJV matrix, which can be efficiently computed (the 
number of operations required is polynomial in the size of 
the chain). This formula is just a different version of the 
one used in [6[, where the Loschmidt echo was written 
in terms of the two-point correlators of the environment 
chain; our formula is improved as the matrix in [6J] was 
2N x 2N. 

In what follows we will analyze the decoherence pro- 
cess for the two-qubit system. For this purpose we would 
need to compute L a b. c d for all values of the labels (ab) 
and (cd) . The formula we presented above enables us to 
compute this echo when either (ab) or (cd) are equal to 
(00) (provided that the initial state of the environment is 
an eigcnstatc of the Hamiltonian Hqq). In order to study 
the decay of the other density matrix elements, we need 
to compute the survival probability of the initial state 
\Eq) evolving first with Hamiltonian H a \, and then evolv- 
ing backwards with H cc i- So far, we do not have a simple 
formula like (fl"12|) enabling us to efficiently evaluate the 
echo in the most general case (see discussion in the Ap- 
pendix). Thus, we will restrict ourselves to consider the 
evolution of initial states for which equation (TT2]) provides 
the full answer. These initial states are a superposition 
of two computational states of the system, one of which 
is the state 1 00) . The case where the other state is either 
1 01) or 1 10) corresponds to a situation in which one of the 
qubits remains in state |0), without interacting with the 
environment, while the other is affected by the chain in 
a way that has been extensively treated in [f| . 

In order to study the effects of the environment over 
the entanglement in the system, we shall take the initial 
state to be of the form \$) = a\00) + (3\11). It should 
be noted, however, that this case is also reducible to a 
single-qubit system, interacting with two separate sites of 
the chain (which has been partially studied in Q). This 
equivalence is illustrated in Figure [1] But even if the 
two problems are equivalent, the case of separate qubits 
interacting locally with distant sites is a more realistic 
situation. It also allows us to address new interesting 
questions related to the entanglement decay induced by 
the environment, and its dependence on the distance be- 
tween the sites to which the two qubits are coupled. In 
this way we will also generalize previous studies where 
the two qubits were treated as "central qubits" (each one 
interacting with equal strength with all the sites in the 
chain), as in 0]. 



B. Disentanglement and decoherence 

Both the entanglement within the two-qubit system 
and that of the system with the chain can be analyzed by 
studying the Loschmidt echo (fT^]) . Because of the initial 
state chosen, the only relevant echo we need to compute 
is L = Loo, n- The entanglement generated between the 
system and the chain can be measured from the purity 
of the reduced density matrix of the system: while the 
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FIG. 1: The one-qubit system interacting with distant sites 
of a spin chain is equivalent to the problem of a two-qubit 
system prepared in the symmetric state q|00) where 
the qubits interact locally with distant sites of the chain but 
not with each other. 



initial state of the system is pure, as a consequence of the 
interaction the qubits become entangled with the chain 
and the state of the system becomes mixed. The purity 
as a function of time can be computed as: 

Tr(p 2 (t)) = l-2\a(3\\l-L). (14) 

The interaction with the environment will also reduces 
the entanglement between the two qubits. This can be 
most easily measured by the negativity Af [10( , obtained 
by partially transposing the reduced density matrix of the 
system, and computing the absolute sum of its negative 
eigenvalues (a nonvanishing negativity is an indicator of 
non-separability of the two-qubit state). For the initial 
state proposed, the negativity is: 

N 4 ,(t) = \aP\L 1 ' 2 (15) 

If the initial state contains any entanglement (e.g. 
\a(3\ ^ 0) then the entanglement in the system decays 
as L 1 / 2 , and vanishes only when the off-diagonal terms 
vanish. 

The relation between the decay of the off-diagonal 
terms and the disentanglement of the two qubits is not 
always so straightforward. In fact, if the initial state is 
a mixture of the form p = (1 — p)\tp)(ip\ + p 1/4, the 
negativity is given by: 

Af(t) = max{0, (1 -p)A^(t) -p/4}. (16) 

This means that, for a given value of p, the state becomes 
disentangled when the echo is below a threshold value 
given by L Um (p) = [4|a/3|(-l + l/p)]~ 2 . This is a well- 
known phenomenon that has been called "sudden death 
of entanglement" 0]: the entanglement can vanish at 
finite time even if coherence decays asymptotically (or 
does not decay completely) . This can happen here if we 
mix the initial state with the identity, taking p large 
enough for L to reach the threshold, but small enough to 
have some entanglement at the beginning. 
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C. From mutually indistinguishable to far-apart 
qubits 

In the following we shall analyze the behaviour of the 
echo for different values of the distance d = n— 1 between 
the sites that interact with the qubits. More precisely, 
one of our interests is the study of two distinct limits and 
the transition between them: first, if both qubits are cou- 
pled to the environment through the same site, then they 
cannot be told apart by the environment. We shall call 
this the case of mutually indistinguishable qubits. The 
other limit is the one in which the qubits are coupled 
to the environment through very distant sites. In such 
a case, if the environment is sufficiently large, its effect 
should be equivalent to the one produced by two identi- 
cal independent environments (as reported in Q). Our 
study will enable us to see the transition between these 
two limits. It is worth pointing out that for the indepen- 
dent environment regime to exist the back-action of the 
qubits on the environment must be small. In the case of 
strong back-action, one expects the environment to in- 
duce an effective interaction between the qubits. Thus, 
the long distance limit would be identifiable properly 
in the weak coupling case, whereas the strong coupling 
regime would be plagued with effects associated with 
such environment-induced interactions (which would de- 
pend on the distance between the qubits). 

It is interesting to describe the expected behaviour for 
the echo in the limits mentioned above. Let us consider 
first the case d = where both qubits are coupled to 
the same site. If the coupling strength of each qubit is 
g, then the resulting echo Loo,ii(ffj£) is exactly the one 
corresponding to a single qubit interacting with one site 
of the chain (which we denote as Lo,i)i with twice the 
interaction strength: 



d = 



£oo,ii(ff>*) = L ,i(2g,t). 



(17) 



The other limit is also easy to understand. If the two 
qubits are coupled to distant sites of a sufficiently large 
environment, we expect that each qubit will effectively 
be interacting with an independent environment (as men- 
tioned above, this is true provided that the interaction 
between qubits mediated by the environment is indeed 
negligible). For independent environments it is easy to 
compute the echo for the two qubits (Loo.n (<?, t)) in terms 
of the single-qubit echo (L 0j i(g,i)), as: 



d» 1 



Loo,n(g,t) -> Ll Jg,t). 



(18) 



We will investigate the decoherence in the above lim- 
its and the transition between both, analyzing how the 
entanglement between the qubits is affected by the fact 
that they are interacting with the same environment. We 
will show that the effective interaction induced between 
the qubits in the strong coupling regime manifests in the 
form of revivals of the echo that depend on the distance 
between the interaction sites. 



III. WEAK COUPLING 

Let us consider first the weak coupling case (i.e., 
g -C 1). We will study the echo for different chain Hamil- 
tonians and distances between the perturbed sites, con- 
sidering a representative case g = 0.1. Figure shows 
the echo as a function of time for 7 = 1 and for several 
values of A and d. 
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FIG. 2: The echo as a function of time for two qubits weakly 
interacting with two different sites of a chain with N = 100, 
7 = 1, g = 0.1. The plots correspond, from top to bottom, 
to the cases A = 0, 0.5, 1, and 1.5. The distance between 
sites is d = (full), 1 (dashed), 2 (dotted), and 3 (dash- 
dotted). In the critical case A = 1 we include also d = 4 

(• ■ — ), and 10 ( •); these curves are not shown in the 

other plots because they are intertwined with the others. For 
comparison we include in grey the limit of two independent 
environments (each of the qubits interacting with a separate 
chain of length N). 



The case 7 = 1, A = is peculiar as it can be solved 
analytically. The ground state for the chain corresponds 
simply to an eigenstate of all the Pauli operators Xj , and 
is therefore not entangled. The resulting echo is the same 



for all distances d > 1. The behaviour for A ^ is differ- 
ent depending on whether A < 1, A = 1, or A > 1. In the 
first case the echo oscillates with a frequency of order 1 
and small amplitude; the oscillation is about a value that 
does not change significantly as d is varied (for d > 1). 
In the critical case the echo decays rapidly up to times of 
order unity and logarithmically afterwards, while oscil- 
lating with very small amplitude and a frequency of the 
same order as before. For a given fixed time, the value of 
the echo increases with d, which means the further the in- 
teraction sites, the weaker the decoherence caused by the 
environment. For A > 1 (but still of order 1) the decoher- 
ence, quantified by T> = 1 — L, is an order of magnitude 
smaller, and the initial oscillations decay rapidly so that 
the echo reaches, for times of order 1, an almost constant 
value. This value increases with distance, but unlike the 
critical case, it rapidly saturates. Thus, we observe that 
for 7 = 1, the long distance regime is reached at a rela- 
tively small number of sites: d > 4 seems to be sufficient 
for the non-critical cases. It is worth mentioning that we 
are analyzing the echo for times shorter than the ones 
where finite size effects become important; these effects, 
that appear in the form of revivals, will be considered 
later. 

A similar study is presented for the case 7 = 0.1 in 
Figure [31 This case is much less regular than 7=1. 
For A < 1 the echo oscillates with a frequency of order 1 
and another one an order of magnitude smaller (but still 
different from g); the behaviour with distance is very 
irregular. In the critical case, the echo decays almost 
linearly for times up to t ~ 10 (after which a logarithmic 
decay starts); at short times decoherence decreases with 
distance, but at longer times it is not monotonic: as d is 
varied between and N/2 with t fixed, V increases with 
d for small distances, and then it starts decreasing again. 
For A > 1, we recover a behaviour similar to the case 
7 = 1, but T> is two orders of magnitude smaller (for the 
other values of A the order of magnitude of decoherence 
was similar for 7 = 1 and 7 = 0.1). The long distance 
limit is not reached as rapidly as for 7 = 1. 

In almost all the non-critical cases analyzed, the inter- 
action with the same site corresponds to the strongest 
decoherence (we found exceptions to this rule in cases 
not displayed, as 7 = 0.1, A = 0.9). The reason for this 
result, that at first sight may appear strange, is that the 
coupling is indeed weak. A naive argument supporting 
this is the following: at short times one expects a Gaus- 
sian decay of the echo, of the form exp(— at 2 ). The time 
width of this decay is inversely proportional to ^fa tx g. 
For d = we have a ~ (2g) 2 while in the case of inde- 
pendent environments a cx 2g 2 . As a consequence, more 
decoherence is expected in the weak coupling regime for 
d = 0. 

The previous study shows that environment-assisted 
interactions between the qubits are not strong. In the 
strong coupling regime, that will be studied below, such 
effects are much more important. However, it is interest- 
ing to see that some behaviour of this type arises in the 
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FIG. 3: The echo as a function of time for two qubits weakly 
interacting with two different sites of a chain with N = 100, 
7 = 0.1, g = 0.1. The plots correspond, from top to bottom, 
to the cases A = 0, 0.5, 1, and 1.5. The distance between sites 
is d = (full), 1 (dashed), 2 (dotted), 3 (dash-dotted), and 
4 (• • — ). In the critical case A = 1 we include also d — 10 

( ■); for the sake of clarity this curve is not shown in the 

other plots. For comparison we include in grey the limit of 
two independent environments (each of the qubits interacting 
with a separate chain of length N). 



weak-coupling scenario when studying the evolution for 
long times. In fact, over those longer times one expects 
the finite size of the environment to play a role and some 
recoherence to occur. The existence of revivals in the 
echo was indeed shown in Q for the case of a single spin 
coupled to one site of an Ising chain. The revival time 
depends on the chain parameters: in the case discussed 
in [6[ the revival was found at t ~ N/2 for 7 = 1, A > 1, 
and at longer times for A = 0.9. In the two-qubit case 
we are analyzing here, as the perturbation acts on two 
different sites of the chain, the revivals are related not 
only to the finite size of the chain but also to the finite 
distance between these sites. This can be seen in Figure 
ID there we compare the case d = N/2 (interaction sites 
that are opposite in the chain) with the situation where 



the two qubits interact each with one site of two different 
chains. 




FIG. 4: The echo for two qubits weakly interacting with two 
different sites in the same chain at long times has revivals 
due to finite-size effects. We show the echo for TV = 100 
and distance d = N/2 (full), and compare it with the case 
where each qubit interacts with one site in two separate chains 
of length N (dashed) and N/2 (dotted). The parameters of 
the chain Hamiltonian correspond to the critical case A = 1, 
7 = 1, and g = 0.1. 



Figure [4] is interesting since it shows that some of the 
revival peaks are clearly induced by the existence of the 
second interaction site (and not by the finite size of the 
chain itself). However, the large peak amplitude is due 
to the fact that we are looking at the very special case 
d = N/2. If we decrease the distance, the first peak be- 
comes fainter, while the second one remains the same 
(Figure [5]). Thus, the peak that is associated to the dis- 
tance between the qubits tends to disappear when this 
distance is much smaller than the size of the environ- 
ment (Figure [5] suggests that d < N/A is enough). We 
can then picture the qubits as evolving almost indepen- 
dently, even for long times. This will not be true in the 
regime of strong perturbations, which will be treated in 
the next section. 



IV. STRONG COUPLING 

We consider now the strong perturbation regime, cor- 
responding to a coupling g_^> 1 (we shall take g = 50). 
This limit was studied in [ll[ in connection with the 
problem of universality, namely, the fact that for large 
enough g the echo displays a fast oscillation (with fre- 
quency of order g) with an envelope which is independent 
of g. 

In the regime of strong coupling we expect for qubits 
interacting with the same site a weaker decoherence than 




FIG. 5: The echo as a function of time for two qubits weakly 
interacting with two different sites in the same chain, with 
N = 100, at distance d = 50 (full), 45 (dashed), 40 (dash- 
dotted) and 30 (dotted). The parameters of the chain Hamil- 
tonian correspond to the critical case A = 1, 7 = 1, and 
g = 0.1. The second peak is a consequence of the finite size 
of the chain. The first peak is associated to the distance be- 
tween the qubits: it reaches its maximum for d = N/2, and 
becomes negligible for small enough d. 



for two very distant qubits (behaving as if they were 
coupled to independent environments). This can be ex- 
plained by the universal behaviour we just mentioned: in 
the case d = we obtain Loo,n(<?j t) = Lo t i(2g,t), but 
this is of the same order as Lo j i(g,t) because for large 
enough g the envelope is independent of g. On the other 
hand, in the long distance limit £00,11 (ff, — L% i(g,t), 
and this is smaller than Lo^(g,t) because the echo is 
smaller than 1. The comparison between d ^> 1 and 
d = 0, illustrated in Figure [5] for 7 = 1, A = 0.99, thus 
leads to a result which is quite the opposite of the one 
obtained for weak coupling. 

There is in this regime another interesting effect, as- 
sociated with the fact that the two qubits interact with 
the same chain. In Figure [6] we see that, apart from the 
fast oscillation and the expected decay of the envelope 
(that are already present in the single qubit case), for 
d = 2 there is a beating, or, more appropriately, an echo 
revival. This is a consequence of the presence of two 
interaction sites; indeed, the timescale of the revival is 
related to the distance between these two sites. This ef- 
fect is not important in the weak coupling limit studied 
before, and it is not related to the finite size of the chain; 
it can instead be interpreted as an interaction between 
the qubits through the modes of the perturbed chain. 

The revivals can be analyzed by looking at the spec- 
trum of the Hamiltonian contained in the matrix A, and 
the mixing of creation and annihilation operators given 
by the matrix h, which appear in the main formula (|12p . 
The unperturbed Hamiltonian has excitations with ener- 
gies given by E k = 2[1 + A 2 + 2A cos(27rfc//V)] 1 / 2 , which lie 
between 2|1 — A| and 2|1 + A|. There is a two-fold degen- 
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FIG. 6: Envelope of the echo as a function of time for two 
qubits strongly interacting with two different sites of a chain 
with N = 100, A = 0.99, 7 = 1, g = 50. The plots correspond 
to the cases of distance d = (full black), 2 (dotted), and 
long-distance limit (dashed). For the case d = a part of 
the fast oscillation is included. The shape of the envelope is 
determined by the chain dynamics and the distance between 
the interaction sites. Here, in contrast to the weak coupling 
scenario, decoherence is stronger in the long distance limit, as 
the envelope for large d corresponds to the square of the one 
for d = 0. 



quasiperiodic behaviour of the echo, with falls and re- 
vivals of the envelope governed by the modes of the chain 
contained between the qubits. As mentioned in Section 
HT1 for initial mixed states the echo decay can lead to sud- 
den death of the entanglement. In the same way, now the 
beating associated to the interaction through the chain 
can also lead to a sudden revival of this entanglement. 

Let us now consider the time of the revival, t r , and the 
amplitude of the peak, L r , as functions of the distance d 
between the interaction sites. We consider first the case 
A = 0.99, j = 1, N = 100, which is close enough to the 
phase transition for i?oo and H\\ to belong to different 
phases. The echo for d — has no revivals until effects 
due to the finite size of the chain appear. For d > 1, the 
revival time increases with distance, while the height of 
the peak decreases. In each case, new peaks appear at 
multiples of t r . Figure [7] displays the dependence on dis- 
tance of time and height of the first revival: the revival 
time is linear, which suggests a constant-velocity propa- 
gation of signals in the chain, while the value of L r decays 
as a power law, L r oc d' 1 / 4 . The Fi gure exhibits also a 
peculiar feature associated with the parity of d: when d is 
even the echo is slightly increased over the general trend, 
and it gets smaller for odd d. For distances d > N/5 the 
effects associated with the finite size of the chain appear 
altering the regular trend. 



eracy associated to ±fc symmetry in all but the lowest- 
and highest-energy excitations. These particles are de- 
localized in terms of the original fermion operators c, 
(which are themselves non-local with respect to the spin 
operators, but whose population is given by the one-site 
magnetization ifTU)) ). When the strong effective external 
field is applied in two sites, two eigenvalues of order g ap- 
pear, which are responsible for the fast oscillation of the 
echo. These excitations are associated to combinations 
of the original fermion operators in those two sites. The 
remaining eigenvalues are still of the same order as before 
(but the degeneracy is broken). These other excitations 
are wave-like, but not entirely delocalized: they can be 
split in two groups, corresponding to excitations between 
the interaction sites, and outside this region. The matrix 
h is related to the populations of the particles that diag- 
onalize the perturbed Hamiltonian in the ground state of 
the unperturbed chain. These can, using equation (fT5)) . 
be written as: 

(EoW J ill) rii 11) \Eo) = (h t h) ii . (19) 

We find that the most populated levels correspond to the 
lowest energy excitation (occupying the outside region), 
the excitations in the interaction sites, and those lying in 
the inside region. The "localized" excitations are associ- 
ated to the rapid echo oscillation, and the lowest energy 
excitation to an oscillation with time scale growing as N, 
while the beating in the echo is given by the frequency 
of the lowest energy excitation in the inside region. 
If we take N — > 00 but keep d fixed, we will find a 



60 




FIG. 7: Time of revival and value of the peak as a function of 
distance, for an Ising chain with N = 100, A = 0.99, g = 50. 
The revivals are determined by the chain dynamics and the 
distance between the interaction sites; the revival time is lin- 
ear with d with slope approximately equal to 2, while the 
height of the revival decreases as gT 1 / 4 ; at d > N/5, finite- 
size effects appear inducing alterations of this behaviour. 

The revival peaks are not a consequence of the phase 
transition in the environment; they also appear for values 
of A < 1 such that no phase transition is involved. But 
for A > 1 the revival is no longer there, as a consequence 
of the change in the ground state properties. The times 
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of revival and the values of the peaks are shown in Figure 
[5] for A between 0.5 and 0.999, with d — 3 (maintaining 
7= 1, JV= 100). 



the peaks are less well-dehned for large d and small 7, so 
that these values are merely estimates. 
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FIG. 8: Time of revival and value of the peak as a function 
of lambda, for N = 100, d = 3, g = 50. A is varied between 
0.5 and 0.999; for A > 1 the revival disappears. 

The power-law decay of the revival peaks in the case 
A = 0.99 might seem to be related to the way in which 
the correlations decay in the spin chain. In fact, for the 
critical Ising chain one can show that the slowest decay, 
corresponding to (XiX i+ d), goes as d~ 1 ^ [12[ lilj]. But 
examining the echo far away from the phase transition 
(we took A = 0.5) we found a slower decay of the revival 
peak with d (even when the correlations decay faster). 
Furthermore, the dependence of t r with d in this case is 
not linear, but grows exponentially. Thus, t r seems to 
be related to transport properties of the chain that are 
modified by varying the external field. This should serve 
as a warning not to take too seriously the simple image of 
revivals as the manifestation of a spin wave propagating 
with constant velocity along the chain (the energies of the 
modes of the chain between the qubits do not generally 
scale as 1/d). 

We might also consider changes in the anisotropy pa- 
rameter 7; as an example, in Figure [5] we compare 
the envelopes for 7 = 1 and 7 = 0.1 (with d = 3, 
g = 50, A = 0.99, and N = 100). The shapes of the 
envelopes are different, and also the heights and times 
of the revival peaks differ. We studied the revival time 
t r and the value of L r as functions of the distance d for 
A = 0.99 and different values of 7 between 0.1 and 1. We 
found that once more, the dependence of t r with d is not 
generally linear; for 7 = 0.1 it becomes approximately 
quadratic. Besides, the height of the revival peaks de- 
creases with 7. The dependence of L r with d maintains 
the power-law decay but with a power that varies from 
approximately -0.25 for 7 = 1 to something between —0.3 
and —0.35 for 7 = 0.1. Anyway, it should be noted that 




FIG. 9: Envelope of the echo as a function of time for two 
qubits strongly interacting with two different sites at distance 
d = 3 of a chain with N = 100, A = 0.99, g = 50. The plots 
correspond to the cases 7 = 0.1 (full) and 1 (dotted). The 
shape of the envelope is determined by the chain dynamics 
and the distance between the interaction sites. 



V. CONCLUSIONS 

Here we will summarize our results. We have stud- 
ied the evolution of a system of two qubits that inter- 
act locally with different sites of an XY chain acting as 
an environment. For this purpose we have analyzed the 
Loschmidt echo, which is the square modulus of the factor 
that determines the reduction of the off-diagonal terms 
in the reduced density matrix of the system. Following 
the ideas in [6| and [14j |. we have derived a formula for 
the evolution of the echo which can be evaluated effi- 
ciently. The entanglement within the system and that of 
the system with the chain were connected to the value of 
this echo. We found that for the initial state proposed, 
= a|00> + /3|11) (with a(3 ± 0), the entanglement 
in the system only vanishes when the off-diagonal terms 
are completely suppressed. This is not the case if we 
mix the initial state \tp) with the identity; then the en- 
tanglement will reach zero once that the echo is below a 
certain value that is related to the amount of mixing and 
the coefficients a, (3. 

It was our aim to study the differences in the echo as 
the distance between the qubits was varied, ranging from 
the case when both interacted with the same site to the 
opposite situation when they were very far apart. The 
cases we have analyzed show a rich variety of results. In 
the weak perturbation regime, we notice a very fast ap- 
proach to the long distance limit, and observe that deco- 
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herence typically decreases with distance. Strong system- 
environment interactions, on the contrary, give rise to a 
different kind of behaviour, in which decoherence is larger 
for long distance. This strong coupling regime is charac- 
terized by fast oscillations with a short-time decay quite 
independent of the distance between qubits, but with a 
beating that provokes quasi-periodic revivals of the en- 
tanglement in the system. This effect can be interpreted 
as an interaction of the qubits through the chain, and 
its time scale is determined by the modes of the chain in 
the region between the qubits. The dependence of the 
revival time with distance is different according to the 
values of the parameters in the chain Hamiltonian: we 
have observed a linear dependence for the almost-critical 
case 7 = 1,A = 0.99, but also quadratic and exponen- 
tial growths for other values of 7 and A. The beatings 
we found correspond to A < 1; when the value of A was 
increased over the critical value (with 7=1 fixed), the 
observed revivals disappeared. 

We note that in the strong coupling regime the fast 
oscillation of the echo continually provokes decays and 
revivals of the entanglement between qubits. This is just 
a dynamical transfer of this entanglement back and forth 
from the system to the environment, with a frequency 
given by the coupling constant g. The true loss of en- 
tanglement is produced by the decay of the oscillation. 
If the initial state is mixed, this decay can lead to true 
sudden death of entanglement. When the distance be- 
tween qubits is short, the sudden death can be followed 
by a sudden revival, as the chain allows for mediated in- 
teractions between the qubits. For very long distances, 
entanglement loss becomes irreversible, unless we con- 
sider times which are long enough for the finite size of 
the environment to become manifest. 



VI. APPENDIX 

In this section we show how the formula (TT2"|) for the 
echo can be obtained. Taking into account the relation 
(fT3")) between the operators that diagonalize the Hamil- 
tonians, 



V 



(0) 



k 



the two different vacuum states \Eq) (for the unperturbed 
Hamiltonian) and \E' ) (for the perturbed one) can be 
connected by [la] : 



\K) 



(20) 



with G = — g _1 h (for the sake of simplicity, we shall first 
assume that g is invertible, and sketch the most general 
case afterwards). The echo can then be calculated from: 

L(t) = \(E \e- iH ^\E )\ 2 

cx {{E'ole-^e-^^e^^lE'v)] 2 (21) 



where in the last expression A is a diagonal matrix con- 
taining the energies corresponding to the different parti- 
cles, and all super-indices are omitted as all states and 
matrices refer to the perturbed Hamiltonian H\. By 
introducing two identities in terms of fermionic coher- 
ent states between the exponentials and integrating two 
times, we obtain Hal: 



L(t) cx 



dai . . . da.Nd(3i . . . d/3 



N 



cxp { ^dGd + (3*e ltA d - ~/3*G/3* 



(22) 



where a, f3* are Grassman iV-tuples. This is a Gaussian 
integral, that can be solved to: 



L(t) cx 



det 



At A 



-e itA 
-G 



(23) 



Using properties of the determinant and the facts that A 
is diagonal and g ± h is orthogonal, and imposing that 
L(0) = 1, the result can be rewritten in its final form 
(HI: 



L(t) = |det(g + he iAt )| 

In case g is not invertible, the relation ((20)) between 
the vacuum states can be generalized by using inter- 
mediate sets of operators in the following way: by the 
singular value decomposition, g = UDV with D diag- 
onal, U, V orthogonal. We define new fermionic opera- 
tors ef = E k U kj V { k °\ ^ = Y, k VjkV k l) - The linear 
transformation between the £M, £^0) now nas ]j instead 
of g, and the vacuum states are the same because the 
transformation does not mix creation and annihilation 
operators. We can assume Dj = for j < ji, and the 
remaining eigenvalues to be nonzero. By interchanging 

particles with antiparticles (£^°' ) <-* ^j' ') for every index 
such that Dj = we obtain a linear transformation with 
an invertible matrix. The calculation of the echo follows 
the same steps as before, except that it is necessary to 
treat indices j < j\ separately. After the Gaussian inte- 
gration, we are left with an expression of the form (|23p 
in which only indices j > ji appear. The desired result 
can be achieved by conveniently introducing some rows 
and columns in the matrix, in such a way to include the 
parts of the matrices with 7 < j\ without changing the 
value of the determinant [14| . 

The present derivation cannot be easily extended to 
more general cases, with two different evolution opera- 
tors as in ([7]), for in that case there are three sets of 
operators that must be related. The steps taken in the 
above calculation then lead to a result that involves a 
sum over an exponential number of terms, and so cannot 
be efficiently evaluated. 
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